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Abstract. We outline an elementary method for proving numerical hypergeo- 
metric identities, in particular, Ramanujan-type identities for I/tt. The princi- 

■ pal idea is using algebraic transformations of arithmetic hypergeometric series to 
^ I translate non-singular points into singular ones, where the required constants can 

■ be computed using asymptotic analysis. 

(N , 

... I have realised that a good translator deserves honour in our literary 
^ I circles, because he is not a craftsman, not copyist, but the artist. He does 

■ not photograph the original, as was commonly believed then, but recreates it 

creatively. The text is the original material for his complex and often inspired 
creativity. Translator is first of all talent. To translate Balzac he needs to 

I at least partially transform into Balzac, to assimilate his temper, to get his 

enthusiasm, his poetic sense of life. 

K. I. Chukovsky, a High Art 

(N 

> . 

QQ . 1. Introduction 

x/^ [ In what follows we keep the standard notation 



O 
(N 



T{s + n) 



O , ^ ^" r(.) 

for Pochhammer's symbol (also known as shifted factorial). When is a positive 
integer we have {s)n = Y[k=oi^ + particular, (1)„ = n! . 

^ ' The quantity vr is linked to the mathematical legacy of Srinivasa Ramanujan in 

man ' ti-i ^ . r,^ ,^ r ^^■ .■ r .■ 

jan 
like 



^ ■ many different ways. In his development of the theory of elliptic functions, Ramanu- 

■ - -' jan came up [28] with a bunch of computationally efficient representations of l/vr, 




vr 



n=0 



42n) ^ = — , (2) 

^ 26" TT ^ ^ 
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•n—n ' \ / 



547r 



(3) 



n=0 

^ (l)n(i)n.(!)n ^3 ^ ^g^) { = (4) 



n=0 

°° /'3^ 



n!3 ' V 48; TTv^ 



E 

r!,=0 



n!3 ' 'V 882-/ TT 



.2MV4MV4M(^^23 + 21460n)( ^1 =1:^, (5) 



, 1 9^2 
1 + lOr.) = (6) 



^^13 ' "-^^ 34n 

n=0 

°° /'l^ Ci^ 

n\ 



n=0 



f:iiM%ii)i(3 + 40„)i- = i£^. (7) 

n=0 

2yn(i)n(!)» qio3 + 26390n) 



n=0 

°° I'i^ 1 nn2 



(|)n(|)n(|)n,,,_ , _ , 1 99^ 



n!3 ' ^994^^ 271^2' 



these are equations (28), (29), (34), (36), (39), (41), (42) and (44) on the hst in [28]. 
Ramanujan did not hide his interest in computing vr; his comment in [28] about 
identity (8) says "The last series (44) is extremely rapidly convergent." 

The identities do not look hard. In spite of this, their first proofs were only ob- 
tained in the 1980s by the Borweins [7] and by the Chudnovskys [12]. A historical 
account of contemporary techniques for proving Ramanujan's and Ramanujan-type 
formulae as well as their supercongruence relatives can be found in [4, 36, 37]. The 
dominating method which, for example, works for any formulae in (l)-(8) is based on 
modular parameterisations of the underlying hypergeometric series. This modular 
technique cannot be counted as elementary, and it does not cover certain impor- 
tant generalisations of Ramanujan's identities such as recently discovered by the 
first author [18, 19, 20] formulae for I/tt^. Guillera's approach to prove the latter 
is by using the powerful Wilf-Zeilberger (WZ) machinery, which itself is quite ele- 
mentary in nature. Finally, changes of variables allow us to "translate" the already 
derived identities into some harder ones [10, 29, 35, 38]; for example. Bailey's cubic 
transformation [17, eq. (5.6)] 

E ^l^" ^ 2 (|)»(|)n(|)n / 27a:^ \ 

n!3 (4_a;)i/2Z^ ^13 l(4_^)3i 

n=0 ^ ' n=0 ^ ^ ' 



with the choice x = 1/2^ translates the WZ provable equation (2) into the WZ 
resistant equation (3). Note, however, that a starting point of any such translation 
has to be an already established one, by modular or WZ means. 

The principal aim of this note is to demonstrate that in many, if not all, cases 
we may start from a "semi- Ramanujan" identity for l/vr (or for a similar quan- 
tity) whose existence only depends on the asymptotic behaviour of a solution of a 
Picard-Fuchs equation at its singular point. Then translation brings the identity 
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into "traditional" Ramanujan-type formulae. Our approach is illustrated by sev- 
eral examples, in particular, we prove equations (1), (3) (hence (2) by translating 
via (9)), (4), (6) and (7), while equations (5) and (8) (of degrees 37 and 29, respec- 
tively) remain beyond the reach of our techniques exposed here. Ramanujan's and 
Ramanujan-type identities for l/vr we discuss in the paper are all classified accord- 
ing to level and degree [9, 16]. We do not explicitly address the two characteristics 
(though the level does appear below) as they are deeply related to modular (hence 
non-elementary!) parameterisations of the formulae. But it is worth mentioning 
that all identities shown in the paper correspond to relatively low degrees. The 
principal difficulty in dealing with high-degree Ramanujan-type identities is lack of 
efficient asymptotic techniques; see Example 6 and our speculative derivation of the 
conjectural equation (29) in Section 5. 

The rest of the paper is organised as follows. In Section 2 we outline the method, 
while Sections 3 and 4 are dedicated to several illustrative examples. Finally, Sec- 
tion 5 highlights further directions of potential applicability of the method, for 
example, for proving some experimentally observed formulae for l/vr^. 



2. Semi-Ramanujan identities 

The classical theorem attributed to Stolz and Cesaro says that if bn is increasing 
and unbounded and 

lim T r- = 7, 

n^oo bn+1 - bn 



then 

On 

7 



lim — 

n— >oo bn 



as well. Take as an example the n-th partial sums a„ = ^^=0 '^^''^'^ analytical 
solution '^nx"' to a linear differential equation, a solution which converges in 

the disc of radius r where r > is a regular singularity of the equation. If Cn ~ 
const ./--"tt,"-! for some a > 0, then a is automatically a local exponent of the 
differential equation at a; = r; assume that the exponent has multiplicity 1. Then 

oo 

c„x" = 7^1 j + higher powers of ^1 j as a; — )■ r~, 

n=0 ^ ^ 

and to find the leading coefficient 7 we can take 6„ = X]fc=o('^)^/^'' n-th partial 
sums of the power series of (1 —x/r)~^ at a: = r, and use the Stolz-Cesaro theorem: 

Eoo r) r> 

^ = n=0(^nX ^ j^f—; (10) 



x/r) " r»-!>oo (a)„/n! ' 



the interchange of the limits is legitimate because the quotient on the left-hand side 
represents a continuous function in the compact interval < x < r. Below we 
typically choose a = ^ in (10), and our starting example, which has motivated the 
method in this paper and whose proof by means of the argument above is due to 
A. Mellit [27], is as follows. 
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Theorem 1. For m a positive integer and Si, . . . ,Sm in the interval < s < 1, the 
following limit is true: 

y r. (^)n(Sl)n(l - Sl)„ ■ ■ ■ - -j^sinVTSj 

lim VI - a; > — n^— n x = \ \ -. (11) 

n=0 j=\ 

Proof. Here we only use the simple fact 

{s)n{l - S)nn _ 1 _ sin TTS 

n^oo n!^ rfs) rfl — s) IT 



which follows from Euler's definition 

T(s) = lim 



n! ^ 



s)n 

of the gamma function. □ 

It is equation (11) which we may call a semi-Ramanujan formula for l/vr. In 
his manuscript [3] A. Aycock, an undergraduate student of Johannes-Gutenberg- 
Universitat in Mainz, uses some known Ramanujan-type series for l/vr and the trans- 
lation method [38] to prove the limit in Theorem 1 when m = 1 and s G {fililil}- 
But because those instances are already known, even in a greater generality, we 
can. . . reverse Aycock's argument! 

Let us do so, using Bailey's cubic transformation [17, eq. (5.3)] 

(1 E#-" = |:^^^#^(-(r^)" (12) 

which is a companion of (9). 

Example 1. Applying 9^ := x— to the both sides of (12) we arrive at 

ax 

n!3 



(1 - 4x)-'/^ Yl (-2^ + (1 - 4x)n)x" 

_l + 8x^ ilUlUl)n ( Tlx \ 



n=0 

where we let x — t- (—1/8)"^. Note that for y = y{x) := —TJxj (1 — 4x)'^ we have 

_ (1 + 8x)2(l -x) 
(1 -4x)3 ' 

therefore 

(l + 8x)(l-a:)V^ ^ (|)n(|)„(|)„ ^Y 27x 



c^(-i/8)+ (l-4x)3/2 ^ n!^ V (l-4x)^ 

r /I \l/2 ^2)"(6)"^6)" n l 

hm 1-y / > =^ 



71=0 
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by Theorem 1. Comparing the hmit with the one for the left-hand side of (13) we 
deduce 



1 n _ t11/2 °° (1)3 



27r 1 - 4x ^ n! 

n=0 

1 ^a)i 1 



-1/8 



n— n ^ ' 



4^2^ n 

* n=0 

a formula for l/vr of Ramanujan type. □ 

Let us try to formalise this approach. Assume we are trying to prove a formula for 
I/tt (or another constant in general) attached to the point y = yo and the generating 
function Y2^=o ^nV^, the latter satisfying an arithmetic differential equation (usually 
of order 3). Then we look for an algebraic substitution y = y{x) and accompanying 
algebraic multiple g = g{x) such that 

oo oo 

5^C„x" = ^7(x)^A„y(x)" (14) 

n=0 n=0 

also satisfies an arithmetic differential equation, and yo = y{r) where r is precisely 
the radius of convergence of ^^q^"^"- Suppose furthermore that we can then 
compute the limit 

lim a/1 - x/rS^ Cnux'' = lim " = 7, 

x^r- ^ n^oo (2)„/n! 

using Theorem 1 or other means. Then a suitable combination Yl'^=o ^n(a + bn)yQ, 
obtained by differentiation and taking the limit as a; ^ r in (14), is an algebraic 
multiple of 7. 

The series we deal with so far are special cases of the hypergeometric series 



m-'^m— 1 



Cll, . . . ) Clm 
bli . . . , bm~i 



^ Sr-^ (aijn • • • [am)n 

^ n\ {bi)n-- ■ {bm.-l)n 



n=0 



which converges in the disc < 1. If all the parameters ai, . . . , and foi, . . . , 6m- 1 
are positive, the function defined by the series can be analytically continued to 
C \ [0, +00) via the Barnes integral. 



m-f^m— 1 



dl , . . . , (3j72 

bl, . . . , bm-1 



X 



T{b,)---Tih 



'm—l ) 



r(ai) ■ ■■T(ar, 



2mJ_^_i^ T{bi+t}---T{bm-i+t} 

where the line of integration Ret = —a separates the poles t = 0, 1, 2, . . . of T{—t) 
from those of T{ai + t) ■ ■ ■ T{am + t); for example, we can take a to be any in the 
range < a < min{ai, . . . , am}- With the hypergeometric function defined by (15) 
in mind we have another (well-known!) family of semi- Ramanujan identities. 
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Theorem 2. Let the parameters ai, . . . , and 61, ... , fcrn.-i of the hypergeometric 
function be positive and am < min{ai, . . . , Om-i} the smallest one in the first group. 
Then 



lim i-xY" 



m—l 



ai, 



X 



m—1 

n 



T{aj - am) T{bj) 



r(aj) T{bj-am)' 



Proof. By the residue theorem, the hmit can be interpreted as the residue of 



r(6i) ■ ■ ■ T{bm-i) r(ai + 1) ■ ■ ■ r(a^ + 1) r(-t) 



— X 



r(ai) ■ ■ ■ T{am) T{bi + t) ■ ■ ■ r(6,„_i + t) 
at its simple pole t = —a^ (cf. [2, eqs. (2.3.12) and (2.4.6)] for the m = 2 case) 
Example 2. Write the transformation (12) in the hypergeometric form 

27x 



(1-4X)^/2^F2(2'^2'^2 

By Theorem 2 we have 



1 1 1 



X 



115 

„_P„( 2' 6' 6 
3^21 ^ ^ 



(1 - Axf 



□ 



(16) 



hm {-yY'%Fj\^\^ y\ 



20F 



V3r(|)3 



and 



lim (-y)V2(i+6^^)3F2f^' ^' i = lim {-yf/\FJ~^\\^ y\ 



2^3 



TT 



where 6y := Z/^- By making the substitution y = —27x/ (l — Ax)^, translating them 
to the left-hand side of (16) and taking the limit x — (1/4) ~ we obtain 



2^ 



n=0 



nP 4" 



3 r(|)3 



and Ramanujan's identity (1), respectively. 



(17) 

□ 



The algebraic transformations we use below are all of modular origin. However, 
as soon as they are given one can verify them directly without any reference to the 
underlying modular functions. This makes our strategy "elementary" in the sense 
"modular free." 



3. Examples of Ramanujan's formulae for I/tt 



Example 3. To prove Ramanujan's identities (4) and (7) we use the transforma- 
tion [29, eq. (3.7)] 



3-^2 



113 
2' 4' 4 
1, 1 



256A; 



9(1 + 3A;)^ 



1 + 3A; 
1 + A;/3 



3-^2 



113 
2' 4' 4 
1, 1 



256A;3 



9(3 + A;)^ 



(18) 
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Apply to the both sides of (18) and pass to the hmit k — )■ (l/9)~ using Theo- 

rem 1 with x = 256A;/(9(1 + 3A;)^) on the left-hand side; the result is equation (7). 
Similarly, equation (4) follows from application of Theorem 2 to the function 



(1+4^^)3^2 2 '^4'^ 



113 
4 



X 



15 3 
p t 2' 4' 4 



3-^2 



1, 1 



X 



(19) 



and then letting k — )■ (—1/3)+. We also have a side identity like in Example 2 
from application of Theorem 2 directly to the 3F2 hypergeometric function on the 
left-hand side of (18); namely, 



°° (-) (^) 



E 

n=0 



nr 



1 

'48 



23/V 

35/41(1)4" 



In a similar way, from the limiting cases z — j- (1/2) and z 
Kummer-Goursat transformation (the square of [30, eq. (25)]) 



(9/8)" of the 



112 

3i^2( 



4z(l - z] 



V9 



'115 

3i^2( 



QAz^{l - z) 
(9-8z)3 



(20) 



we produce the Ramanujan-type formulae 

00 



00 fl\ fl\ f5\ 



Un] 



n=0 



W. 



125 



5^5 
27rv^ 



and 



respectively. 



°° (-) (-) 

M,h[,)n ^^^^^^ 

nr 

n=0 



9 

16 



□ 



Example 4. In many cases, the leading term in the asymptotics of a binomial sum 
expression can be obtained by using Stirling's formula and other standard tools of 
asymptotic analysis. For example, the growth of the sequence Cn = ^"=0 (") is 
determined from [26, Theorem 1]: Cn ~ 2'^"+^/^(7™)~'^/^ as n — 00. On using the 
Stolz-Cesaro theorem we obtain 



lim Vl — 16x CnUx"^ = li 

(1/16)- ^ n- 



lim 



a;-i.(l/16) 



n=0 



C„n/16" _ V2 
{h)n/n\ vr ■ 



(21) 



On the other hand, the generating function Yl'^=o CnX"' can be written hypergeomet- 
rically; namely, applying [24, formula (2-53)], Clausen's formula and [14, Theorem 
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3.6] (see also [15]) we have 



n=0 



1 + fc2) 



4 



1 + 22A; - 6F _ 22k^ + 

■1 13 



X 3F2 



2' 4' 4 
1, 1 



256A;(1 - A;2)(l + A; - - 4A; - 

(1 + + 22fc - 6P - 22P + A;4)4 ) ' ^^^^ 



Apply /c— — to the both sides of (22) and compute the limit as A; — )■ {^/lO~^^A^/E—2 — 

Cl/b 

y/5)~ which corresponds to a; — (1/16)~ in (21) on the left-hand side of the resulting 
identity. This produces the first "elementary" proof of Ramanujan's identity (6), as 
no translation of a WZ provable identity is known for it. 
Another limit as 



k^ 

corresponds to taking 



11 + 5\/5 - V250 + 110V5\ 
2 ) 



_ 256A;(1 - A;2)(l + k- k^f{l -4k- 
^ (1 + P)4(l + 22A; - 6P - 22A;3 + k^y ^ 

Therefore, combining Theorem 2 applied to the function (19) with the transforma- 
tion (22) we deduce the Ramanujan-type formula 

EEC ^'M-h) -h 

n=0 j=0 ^ ^ ^ 



given previously in [15, Theorem 5.2]. Applying Theorem 2 directly to the 3F2 series 
on the right-hand side of (22) we obtain the side evaluation 



EE 

n=0 j=0 



\ 4 / 1 \ 

n\ I 1 \ IT 



J J \ 20J 5i/4r(|)4 

which seems to be new. □ 
Example 5. This time we deal with the level 7 sequence [16] 



j=0 

so that 



hm V 1 - 27x > Cnux = lim ' . = 

x->(l/27)- ^ n^oo i 471 



(23) 
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by the Stolz-Cesaro theorem. It follows from [16, Theorem 3.1] and [8, Lemmas 4.1 
and 4.3] that 



CXJ y 

n=0 ^ 



Vl + 13/1 + 49/^2 



13/i + 49/i2; Vl + 5/i + /i2 



^ n!3 V(l + + + 13/i + 49/i2) i ' ^ ^ 



n=0 



Apply to the both sides of this transformation and compute the limit as /i — )■ 
ah 

(l/7)~ which corresponds to x — (l/27)~ in (23) on the left-hand side. The result 
is Ramanujan's series (3) which, in turn, is equivalent to (2). □ 

Example 6. Let Xl^o^"'-^"' 'with = 1, be the holomorphic solution of the 
differential equation Vy = 0, where 

p := ^3 _ ^(^29 + l){lie^ + 11^ + 5) + x'^ie + 1)(121^2 ^ 242^ + 141) 

-98x3(^ + l)(^ + 2)(2^ + 3) 

and 6 = 6x = x—. The equation is equivalent to a difference equation with poly- 

ax 

nomial in n coefficients, for the sequence itself, and the Birkhoff-Trjitzinsky 
method [32, 34] allows us to determine asymptotic behaviour of the latter: 

. C{1 + 2v/2)2" / Ci C2 , \ 

\ n J 
with the numerical value C = 0.44846956652386 . . . which can be recognised as 

^_ ^184+11^2 

4^2713/2 

(An efficient strategy to compute the number is outlined in [33, Section 3].) 

As D. Zeilberger writes in [34], "The Birkhoff-Trjitzinsky method suffers from 
one drawback. It only does the asymptotics up to a multiplicative constant C. But 
nowadays this is hardly a problem. Just crank-out the first ten thousand terms of 
the sequence using the very recurrence whose asymptotics you are trying to find, 
not forgetting to furnish the few necessary initial conditions, and then estimate 
the constant empirically. If you are lucky, then Maple can recognize it in terms of 
'famous' constants like e and vr, by typing identify (C) ;." 

In order to not only discover the constant C numerically but also to prove it 
rigorously, we need more information about the sequence An, for example, a binomial 
sum expression for its terms. 

The above asymptotics is equivalent to 



lim V(l -4x)(l - 18x + 49x2) VA^nx'' = (25) 



2ti 

n=0 



10 
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where r = 1/ (1+2a/2)^, through the Stolz-Cesaro theorem. The underlying modular 
structure (of level 14) allows us to link the series AnX"' to the (level 7) series 

^n^" from Example 5: 



1 °° / u \" 1 °° / V \' 

;l + ,;)(l + 8^)SJ^"((l+^)(l + 8^;)j =(1^4^ 2 ^4 (1^4^ J 



(26) 



The value r corresponds to the choice V = 2-3/2^ so that equation (25) translates 
into the formula 



n,=0 j=0 



TT 



V7 



of Ramanujan type; it can be also derived from the results in [16, Section 4] (it 

corresponds to taking N = 2, q = e'^'^v^, yj = {1 - l/y/2f and A = (3 - ^2)/? 
in the notation in [16]). With the help of (26), (24) and the known algebraic 
transformations of 3F2 hypergeometric series, like (12) and (20) above, we can obtain 
some more Ramanujan- type formulae for l/vr; for example. 



y (|)n(i)n(|)n _ 3^ ^ 55^) ( L\ 

t^, 'vii + 8v^; 



2^25 + 22^2 



TT 



n=0 

The latter identity can be found in [6, eq. (4.75)]. □ 

Here seems to be a right place to briefly comment about underlying modular 
parameterisations. Any known generating function Yl'^=o ^n^"' which gives rise to 
Ramanujan-type series for l/vr is annihilated by a third order linear differential 
equation V = 9^ + ■ ■ ■ G Q[6'](a;), at least in the genus cases as defined by Conway 
and Norton [13]. Define P{x) G Q{x) to be the formal limit of as 6* — oo (see 
also [1, Section 2]); it is the polynomial 

1 - 2 ■ llx + 121x2 _ 2 . 983,3 = (1 _ - 18x + 49x2) 

in Example 6 above. Then the Conway-Norton modular parameterisation [13, Ta- 
ble 3] x{t) = 1/Jg{t) and the use of the Stolz-Cesaro theorem result in 

00 

hm v/pR E = TT (27) 

n=0 

where i is the level of the corresponding modular group G and r is a suitable (usually, 
the smallest) zero of P{x). This can be counted as an arithmetic generalisation of 
the m = 1 case of Theorem 1; a heuristic analogue of (27) exists for m = 2 as well [1]. 
As hinted in [38], the semi- Ramanujan identity (27) can be further translated (at 
least, in theory) into other Ramanujan-type formulae for l/vr whenever the modular 
parameter r is chosen from the same quadratic field Q(-\/— £). 
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4. Multi-parameter translation 
Example 7. Begin with an identity due to Whipple [17, eq. (2.2)]: 



11 



X 



I + p — q, 1 -\- p — r 



X 



( i + p - g - r, i + ip, I + 
1+p — q, 1+p — r 



3^2 



Ax 



Then apply the operator 6^ and take the limit as x — 
the right-hand side, to obtain the identity 



using Theorem 1 on 



y ^t^^^"^^t!^"^^^'^: (1 + 2p + 4n) (-1) 

n! (1 ' " " ' "~ 



n=0 



(1 + p - g)„(l + p-r)n 
2i-p 



^l)p-(?(l)p-r 



2T{l+p- q)T{l+p-r) 



^ (iWr(|)p/2(t)p/2 r(i + p-g-r)r(i + p) 

This generalisation was first obtained in [11] using a different method. 

In a similar way, the limiting x (—1/8)"^ case of Bailey's identity [17, eq. (5.3)] 



1 + 2p, i - 2g, i + 2q 
p- q, 1+P + q 



X 



1+p - q, 1+P + q 



3-^2 



27x 



(2J 



{1-Ax)\ 

which is a general form of (12), leads us to the following generalised Ramanujan-type 



series: 



^ (l + 2p).(l-2g)„(l + 2g)„. , ^ 



n=0 



+ P - g)n(l + P + g)n(l)n 



25/2+2p^i/2 Y{l+p- q)T{l+p + q) 



r(| + fp)r(| + |p)r(| + |p)' 

which was first obtained in [31] in a different way. 

If instead we compute the limit as x — ?■ (1/4)+ of (28) using Theorem 2 on the 
right-hand side, we find that 



^ (^ + 2p)^(i-2g)„(i + 2g)„ 1 
^ n\{l+p-q)n{l+p + q)n 4^^ 

24(l+p)/3^ 



T{l+p- q)T{l+p + q) 



31+2P r(i + hni + hni + y- qm + b+qy 



which generalises equation (17). 



□ 



Due to the presence of extra parameters, the identities in Example 7 can be 
directly shown by the Wilf-Zeilberger method. 
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Example 8. Take 



and consider the related transformation [35, Theorem 2] 

Ax 



n=0 ^ ^ ^ ^ n=0 



15 



6'a;-Differentiating and taking the limit as x — (— l)"*", we recover the identity 

which is due to Ramanujan. The general case of the transformation [35, Theorem 1] 
results in an identity also recorded by Ramanujan [5, Entry 31, p. 41] which, in 
turn, is a special case of Whipple's transformation. □ 

5. Conclusion 

It seems plausible that a suitable generalisation of the translation techniques 
considered in the paper, maybe an extension of a purely computational nature, will 
make it possible to prove some experimentally observed but not yet proven formulae, 
hke 

g (|)n(|)n(|)n(|)n(i)n ^ ^^^^ ^ ^g^Qn^) i- = (29) 

n=0 



y (l)n(i)n(|)n(|)n(i)n. , V = i| 

n=0 ^ ^ 

°° (1\7 



V ^ (1 + 14n + 76n2 + 168n=^) = ^. 

n=0 

Note the resemblance of these equations with (7), (4) and (2), respectively, and the 
fact that for each of the latter we already know "suitable" transformations from 
Examples 3 and 5. For instance, the proof of (29) might go as follows. 
Define the sequence C„ by the equation 



^ / 256A: 
^ ^9(1 + 3. 



3k)\ 



k/3 



Then the generating function Y2'^=o ^nX^ satisfies a linear differential equation, 
equivalently, the sequence C„ satisfies a difference equation with polynomial co- 
efficients. Assuming we can guess from the equation a hypergeometric (binomial) 
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expression for the terms C„, we can find the asymptotics 

Cn ~ 27671-^/2^-3/2 as n ^ cx), 

so that 

hm \/\ — X CrJix^ = — ^ 

71=0 

by the Stolz-Cesaro theorem. Applying k— to the both sides of (30) and computing 

dk 

the hmit as — )■ (1/9)- give us the desired identity (29). 

Another aspect of the method is its potential applicability to proving so-called 
Ramanujan-type supercongruences [21, 22, 37]. This belief is supported by some 
recent results of L. Long and her collaborators [23, 25]. 
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